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Abstract 

We derive rigorous bounds on corrections to Einstein gravity using unitarity and analyt¬ 
icity of graviton scattering amplitudes. In D > 4 spacetime dimensions, these consistency 
conditions mandate positive coefficients for certain quartic curvature operators. We sys¬ 
tematically enumerate all such positivity bounds in D = 4 and D = 5 before extending to 
D > 6. Afterwards, we derive positivity bounds for supersymmetric operators and verify 
that all of our constraints are satisfied by weakly-coupled string theories. Among quadratic 
curvature operators, we find that the Gauss-Bonnet term in D > 5 is inconsistent unless 
new degrees of freedom enter at the natural cutoff scale dehned by the effective theory. 
Our bounds apply to perturbative ultraviolet completions of gravity. 
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1 Introduction 


Low-energy effective theory describes quanta interacting indirectly through kinematically inac¬ 
cessible states. The dynamics are characterized by an effective action that typically includes all 
interactions permitted by symmetry with couplings of order unity. However, in certain cases, 
self-consistency at long distances imposes non-trivial constraints on the coefficients of effective 
operators. This is famously true in the theory of pions, where the operator coefficients of certain 
higher-derivative operators are required to be strictly positive to ensure causal particle propaga¬ 
tion together with unitarity and analyticity of scattering amplitudes at complex momenta [1-5]. 

As low-energy criteria, causality, unitarity, and analyticity impose constraints that are inde¬ 
pendent of the detailed ultraviolet dynamics. Consequently, these consistency conditions offer 
special utility in the context of quantum gravity, where the ultraviolet completion is not known 
with certainty. For instance, such bounds have been derived for the effective theory of gravitons 
and photons [6], where consistency necessitates large charge-to-mass ratios precisely of the form 
of the weak gravity conjecture [7]. 

Notably, a proper diagnosis of causality violation in curved spacetime is subtle since particle 
propagation can be locally superluminal even in healthy theories. For example, it has long been 
known that photons with certain polarizations travel superluminally in the vicinity of a black 
hole in the effective theory of photons and gravitons describing our actual universe [8]. Instead, 
a more global measure of causality, e.g., the existence of closed timelike curves, is necessary to 
establish a true pathology. On the other hand, unitarity and analyticity offer alternative criteria 
that are mathematically rigorous and applicable in the flat-space limit. 

In this paper, we systematically derive new constraints on curvature corrections in gravity 
from unitarity and analyticity. The graviton effective theory is described by the action^ 

/ OO 

( 1 ) 

n=l 

where £„ are contributions to the action entering at order 2n in the derivative expansion and 
£1 = ^ and C 2 = + R^) (2) 

are the Einstein-Hilbert and Gauss-Bonnet terms. We assume the Gauss-Bonnet form for C 2 
throughout since this is the unique ghost-free quadratic curvature invariant [9,10] in D dimen- 

^We work in mostly -I- signature and write k = VSttG, ^pi,, and R^„pa = 
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( 3 ) 


sions. Moreover, we restrict our analysis to effective theories in which £4 takes the form 

7 

£4 = CiOi, 
i=l 

expressed in terms of the minimal basis of quartic Riemann operators^ in Ref. [12], 

/O _ I? DajSyS p _ Tifiupa p S pa/37 p 

^ -tioif^ryg L/2 .K -^pup ^ (j^Oi0'yS 

Us — K ttpo.75 Ui — n rtpj, /tp^, r£o-/3..y5 

/o _ JD 01/3 T3 ■y S TD pr\ _ rypupa td a j3 t3 1 ^ TD 

Us tx rXp Q, itcryps Uq 11 JrC^ ^ rC^ ^ ^n-yaS 

py _ TDUVpcr TD a 0 TP y S TP 

Uy - ^ p K/s^aS- 

Note that linear dependences arise among operators as the dimension D of spacetime decreases. 
At quadratic order, £2 is unphysical in H < 3, a total derivative in H = 4, and dynamical 
in H > 5. Meanwhile, at quartic order, the number of algebraically independent operators 
Oi in D = 4, 5,6,7, 8 is 2,4, 6 , 6 , 7, respectively [12], with one linear combination—the eight¬ 
dimensional Euler density—a total derivative in 77 = 8 and hence dynamical only in 77 > 9 [13]. 

Our analysis hinges on the on-shell four-point graviton scattering amplitude, M, whose 
forward limit is intimately linked to the total cross-section by well-known analyticity argu¬ 
ments [1,5]. By marginalizing over the external graviton polarizations, we can then systemat¬ 
ically derive a rigorous and inclusive set positivity bounds on the coefficients of operators in 
the graviton effective action. Throughout, we assume a perturbative ultraviolet completion of 
gravity, so there exists a well-dehned expansion in h. 

We begin with an analysis of quartic curvature corrections, proving that in 77 = 4 the 
coefficients of the (77^)^ and (7777)^ operators are positive. Our results precisely match those of 
Ref. [14], which derived bounds from the condition of locally subluminal graviton propagation. 
We then generalize our arguments to 77 = 5 and 77 > 6 . Subsequently, we obtain positivity 
constraints on supersymmetric operators in general 77, which in the literature are sometimes 
denoted tstgR^ and As a consistency check, we verify that all our constraints are 

satisfied in the bosonic, type II, and heterotic string theories. 

Moving on to quadratic curvature corrections, we argue that unitarity and analyticity exclude 
theories in which A S> 1 with no new degrees of freedom present at or below the mass scale 
A ~ the natural cutoff associated with the Gauss-Bonnet term and the derivative 

expansion. Our results precisely accord with those of Maldacena et ah [15], who demonstrated 
that this class of theories is inconsistent with global causality. 

^Applying leading-order equations of motion to £„ is equivalent to a field definition modulo new terms generated 
in £n-i-i. Repeating this procedure at each order, we can freely impose R = Rpv = 0 in a pure gravity theory [11]. 


4 



The remainder of this paper is organized as follows. In Sec. 2, we review the arguments 
of Ref. [1] whereby unitarity and analyticity impose rigorous positivity bounds on operator 
coefficients in effective theories. Afterwards, in Sec. 3 we apply these methods to establish the 
positivity of certain coefficients of quartic curvature operators, starting in D = 4 and D = 5 and 
then generalizing to D > 6. We then apply our bounds to supersymmetric theories and string 
theories. Finally, we study quadratic curvature operators in Sec. 4 and conclude in Sec. 5. 

2 Analyticity Argument 

In this section we review how operator coefficients in effective held theories are constrained 
by the analyticity of scattering amplitudes at complex momenta. Our analysis follows that of 
Ref. [1], which derived positivity bounds on operator coefficients by relating the low-energy limit 
of forward amplitudes to strictly positive integrals of cross-sections. 

Our object of interest is the on-shell amplitude M describing four-point graviton scattering in 
D dimensions. Here, the choice of the external polarizations is built into the functional form of 
M, as is the case for helicity amplitudes in D = 4. From this viewpoint, helicity is just a quantum 
number labeling the external states, no different from baryon or lepton number. Sometimes it 
will be useful to view M as a function of the external particle labels, M = M(l,2,3,4), and 
other times as a function of the kinematic invariants, M = M{s,t,u), where 

s = — [ki k 2 )^ 1 t = —(fci - 1 - ^ 3 )^, u = —[ki ^ 4 )^, (5) 

working in the convention where all momenta are incoming, so /ci -|- /c 2 + + ^4 = 0 . 

To derive constraints from analyticity, we will be interested in scattering amplitudes that 
are simultaneously forward and invariant under crossing in the f-channel. Formally, t-channel 
crossing symmetry implies invariance under swapping particle labels 1 -H- 3 or 2 -H- 4 while 
leaving the functional form of M —which encodes the polarization choices—hxed, so 

M(l, 2, 3,4) = M(3, 2,1,4) = M(l, 4, 3, 2) = M(3,4,1, 2). ( 6 ) 

For external gravitons, crossing symmetry is ensured if the exchanged states are indistinguish¬ 
able. This happens in D = 4 if the states have identical helicity and more generally in D 
dimensions if the states have the same polarization relative to their momenta. Mathematically, 
M is crossing symmetric provided the momentum and polarization of particle 1 are related by 
an improper Lorentz transformation to those of particle 3, and likewise for particles 2 and 4. In 
terms of kinematic invariants, a crossing symmetric amplitude then satishes 

M{s,t) = M{-s (7) 
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where the momenta are swapped but the polarizations relative to momenta—which m. D = A 
are the helicities—are untouched. 

Meanwhile, the forward limit of the amplitude, M{s,t —)■ 0), corresponds to an identihcation 
of particles 1 f-)- 3 and 2 -H- 4. This is achieved simultaneously with crossing symmetry if we 
restrict to following kinematic regime: 


forward and 
crossing symmetric 


-H-(-/ci,ei) and (/c 4 , £ 4 )-H-(-fc 2 , 62 ), 


( 8 ) 


where ei and 62 are real linear polarizations. We choose a basis of linear polarizations because an 
amplitude with hxed external circular polarizations cannot be simultaneously crossing symmetric 
and forward. 

The forward and crossing symmetric amplitude, M{s,t —)■ 0), can then be expanded in a 
power series in s and t. While analytic singularities in s or f arise, their form is severely restricted 
by the locality of the underlying theory. As noted earlier, we assume throughout a perturbative 
ultraviolet completion of gravity, so we are justified in considering only the leading contribution 
in the h expansion, he., tree-level exchange. 

At tree level, analytic singularities in kinematic invariants enter at worst as simple poles. 
Moreover, a t-channel singularity in the forward limit can only arise from non-local terms cor¬ 
responding to graviton exchange induced by the leading Einstein-Hilbert interactions, so the 
general form for the forward amplitude is 

CX) 

M{s,t ^ 0) = fnS'^ + 0{s‘^/t). (9) 

n=0 

The hrst term is regular, as it is generated by heavy particle exchange, while the second term 
is singular because it comes from f-channel graviton exchange scaling as ~ s^/t. The form of 
Eq. (9) together with the crossing symmetry relation of Eq. (7) implies that 


M(s, f —>■ 0 ) = M(—s, t —)■ 0 ) -|- 0{s), 


( 10 ) 


where the hrst term arises because the limit of a regular function is the function evaluated at the 
limit of its arguments, while the second term is a residual contribution from applying crossing 
symmetry to the singular 0{s^/t) contribution. 

The parameters /„ depend on the coefficients of operators in the effective action of the 
low-energy theory. To determine analyticity constraints, we consider the order n residue of 
M(s, t —)■ 0 ) in the complex s plane, yielding 
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where C denotes a small contour encircling the origin. 

As previously noted [1,6], the 0{s^/t) singular contribution is formally inhnite in the strictly 
forward limit and therefore a major obstacle to deriving bounds from analyticity. Nevertheless, 
for n ^ 2 this term is eliminated by the residue theorem. While forward singularities of order 
s^'/t can arise from loop-level graviton exchange diagrams, we are working at leading order in 
the h expansion so these contributions are formally subdominant. On the other hand, n = 2 
is more subtle, but we will show that in certain parameter regimes the 0{s^/t) term can be 
subdominant to the rest of the amplitude, allowing for a bound to be placed. In any case, we 
leave a detailed discussion of these issues for later sections and for now simply drop the 0{s‘^/t) 
contribution. Terms subleading in the forward limit of the Einstein-Hilbert amplitude must by 
power counting go as 0{s) and will thus be eliminated in the contour integral for all n > 2. 

By Cauchy’s theorem, we can blow up C into a new contour C that runs just above and 
below the real s axis, plus a circular boundary contour at infinity, yielding 


fn 


1 

27ri 



ds 

gn+l 


Disc[M(s,t —)■ 0)] -|- boundary integral. 


( 12 ) 


where Sq is any scale above zero and below the hrst massive threshold in the ultraviolet comple¬ 
tion. We note that Disc[M(s,t —)■ 0)] = M{s + ie,t ^ 0) — M{s — ie,t —>■ 0). By the Schwarz 
reflection principle, M{s*,t ^ 0) = M{s,t 0)*, so 


Disc[M(s,t —)■ 0)] = 2ilia[M{s,t —)■ 0)]. (13) 

The crossing symmetry relation in Eq. (10) then implies that 

Disc[M(—s,t 0)] = —Disc[M(s,t —)■ 0)], (14) 

dropping the 0{s) term that arose from the 0{s^/t) singularity. 

Throughout this paper we assume that \M{s)\ is less divergent than Isl"" at large s so that 
the boundary term in Eq. (12) can be dropped.^ This is a physically reasonable assumption 
applicable to any ultraviolet completion in which the large s behavior of the amplitude at fixed 
hnite physical t s grows more slowly in s than the Einstein-Hilbert contribution, which scales 
as s‘^/t. A theory that fails this criterion would actually have worse ultraviolet behavior than 
Einstein gravity. Operationally, this translates into the assumption that \M{s)\ grows more 
slowly than |sp at large s. For example, this can be verihed explicitly in the Regge behavior of 
string theory amplitudes, which scale as where a{t) <2 for t < 0 [1]. 

^Strictly speaking, positivity bounds only require that the boundary term be non-negative, which is sometimes 
true given specific assumptions about the ultraviolet [16]. We do not consider this possibility here. 
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Combining Eq. (12) with Eqs. (13) and (14) yields 


4-1 r°° ds 

fn = ^Im[M(.,t ^ 0)]. (15) 

For n odd, this result trivially implies fn = 0 as required by crossing symmetry of M, while 
for n even, it imposes a positivity condition. In particular, we use the optical theorem to write 
Im[M(s,f —)■ 0)] = S(t(s), where cr(s) is the total cross-section.^ Crucially, in an interacting 
theory with new heavy states, a(s) is strictly positive, so 


fn 


2 

71 



> 0 , 


( 16 ) 


thus establishing positivity of fn for even n. 

Here /„ corresponds to the s" contribution to the low-energy amplitude, which is proportional 
to the operator coefficients of £„. By power counting, we know that the low-energy amplitude 
can be expanded in powers of Mandelstam variables, so 


M = Y,Mn, (17) 

71=1 


where the leading contribution arises from the Einstein-Hilbert action, which in the forward 
limit gives an amplitude 

2 2 

Mi(s,t ^ 0) = ^ 2 pa ^2 —^-1" ^('5); (18) 

where the 0{s) terms are regular in the forward limit. The remaining contributions are 
generated by £„• In the subsequent sections, we derive precise analyticity bounds for the 
quartic and quadratic curvature corrections, £4 and £ 2 - 


3 Bounds on Quartic Curvature Corrections 

In this section we derive bounds on £4, which encodes quartic curvature corrections to Einstein 
gravity. The leading contributions from £4 are quartic graviton vertices, which contribute 
to graviton scattering amplitudes via contact interactions. Since these corrections are free 

^While the total cross-section diverges in the presence of a t-channel singularity, Im M(s, t —>■ 0) and by extension 
cr(s) = Im M(s,t —> 0)/s are really proxies for the finite sum over all residues of heavy states in the complex 
s plane. By factorization, each contribution is positive since M{hh —?> hh) ^ —M(hh —>■ X)M{X —>■ hh)/{s — 
-b id) on the s-channel resonance of a massive state X of mass m. 







from kinematic singularities, their forward limit is regular. Thus, to obtain a forward, crossing 
symmetric amplitude, we simply set t = 0, €3 = ei, and €4 = € 2 , as derived in Eq. ( 8 ). 

After a lengthy calculation, we compute the quartic corrections to the graviton scattering 
amplitude in the forward limit to be 

4 4 

M4(s,t —)■ 0 ) = 2 [(2c6 + ^ 2 pcT^ 2 ^^) + (32ci + 4c2 + 2 c 6 )(ei^,ye 2 ^‘^)^ 

+ ( 4 C 2 + I 6 C 3 + 2 c 6 )(e 4 ^ €-2vp^l^ ^2crp) + ( 4 C 2 + 8 C 4 + 2 C 7 )(e 4 ^ ^lup^ 2 ^ ^‘2.ap)\ ■ 

Eq. (16) bounds / 4 , corresponding to the coefficient of the contribution to the amplitude, to 
be positive. To determine the constraint on the coefficients of £ 4 , we should marginalize over 
all possible values of the independent polarizations, ei and € 2 . 

To determine the full set of bounds, it will be convenient to map the question of positivity to 
a linear algebra problem. To do so, we work in the center-of-mass frame, where the polarization 
tensors, and 62 ^ are real, symmetric {D — 2)-by-(Zl — 2) matrices satisfying the usual 
tracelessness and normalization conditions, 

Tr(ei) = Tr(e 2 ) = 0 and Tr(ei • ei) = Tr(e 2 • € 2 ) = 1. (20) 

Furthermore, we can dehne Hermitian matrices Hj^ = {ei,e 2}/2 and = i[ei,e 2]/2 encoding 
the polarization information, which enter the amplitude in terms of the invariants 

X = Tr(hf+)Tr(i7+), y = Tr(Ef+■ hf+), 2 = Tr(hf_ • hf_). (21) 

We can then express the analyticity bound as 

(2c6+C7) + (32ci+4c2+2c6)x+(8c2 + 16c3 + 8c4 + 2c6 + 2c7)|/+(—I6C3 + 8C4 — 2c6+2c7)z > 0, (22) 

for all {x,y,z) spanned by the graviton polarizations ei and 62 - What is the allowed space of 
{x,y,z)7 An obvious set of necessary conditions are 

0<x,y,z<l and y + z < 1, (23) 

from familiar linear algebra inequalities. In general D, Ending the space spanned by the allowed 
(x, y, z) is a highly nontrivial problem in matrix inequalities. 

In the next subsections, we will study various physically well-motivated scenarios, including 
general theories in H = 4 and supersymmetric theories in arbitrary D. 
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3.1 Theories in = 4 


The number of linearly independent curvature invariants monotonically increases with the di¬ 
mension of spacetime. In D = 4, there are only two independent quartic curvature invariants. 
Hence, £4 in Eq. (4) collapses to 

£4 = ciOi + ciOi, (24) 

where Oi is dehned as in Eq. (4) but Oi is unique to D = 4, 

O, = and O, = R^^^s, (25) 

where R^ypa = RpJ^^eai 3 pa /‘2 is the dual Riemann tensor. The operator Oi can be written as a 
linear combination of any two the operators in Eq. (4) modulo contributions proportional to R 
and Rpui which can be eliminated by the equations of motion. For example, 

Oi = 402 - 4O3 = -4O2 + 8O4 = • • • , ( 26 ) 


corresponding to a choice of operator coefficients, (ci, 4ci, —4ci, 0 , 0 , 0 , 0 ), (ci, —4ci, 0 , 8 ci, 0 , 0 , 0 ), 
etc. The ellipses in Eq. (26) denote equivalent representations in terms of other operators, which 
are not unique due to the linear dependence in O = 4 of all but two of the operators in Eq. (4). 

In O = 4, the invariants (x, y, z) are constructed from real, symmetric, traceless 2-by-2 
matrices, which we can parameterize by 


ei = Cl • (j/\/2 
62 ^ ^2 ■ 


(27) 


where el and el are real unit polarization vectors and a are the Pauli matrices. Since ei and 62 
are real and symmetric, they only have components in ai and a^, since a 2 is imaginary and anti¬ 
symmetric. From standard matrix identities, we see that {ei, e 2 } = £i -£2 and [ei, e 2 ] = i{eixe 2 )-a. 
Dehning 9 to be the angle between ei and £ 2 , we obtain 

{x, y, z) = cos^ 6 (1, |, 0) -h sin^ 9 (0, 0, |), (28) 

which dehnes an interval whose endpoints are (1,4,0) and (0,0,4). Inserting these {x,y,z) 
values, along with the coefficient choice given by Eqs. (24) and (26), the bound (22) takes the 
suggestive form 

Cl cos^ 0-I-Cl sin^ 6^ > 0, (29) 


which obviously implies positivity of both coefficients separately. 


Cl > 0 and ci > 0, 


(30) 


which correspond to parallel or perpendicular polarization vectors, respectively. Our results 
exactly coincide with those derived from requiring subluminal graviton propagation [14]. 
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3.2 Theories in D = 5 


In D = 5, there are four linearly independent quartic curvature invariants. For the sake of 
generality we use the basis of Eq. (4) with the linear dependences among operators assumed. 
For this analysis, we ascertain the physically allowed region for the invariants {x,y,z), which 
in D = 5 are constructed from real, symmetric, traceless 3-by-3 matrices. This requirement 
constrains {x, y, z) to he in the plane 1 + 2x — 6y — 2z = 0. Specihcally, (x, y, z) are restricted 
to a planar triangular region, 

3 

{x,y,z) = ^TiVi, (31) 

i=l 

dehned by three vectors 


ni = (0,0, i), W 2 = (l, |,0), and n 3 = (0,i,0) (32) 

for the real parameters ti,T 2 ,T 3 > 0 such that ti + r 2 + rs = 1. The vertices (32) of this triangle 
can be reached by choices of physical polarizations. In particular. 


Vl 


V 2 


V3 


£l = 


ei = 


ei = 



62 = 


62 = 


62 = 



(33) 


Plugging in Eqs. (31) and (32) back into Eq. (22), we obtain 

Ti(-8C3 + 4c4 + C6 + 2C7) 

+ T2(32ci + 8c2 + 8C3 + 4c 4 + Scg + 2C7) > 0, (34) 

. + ^T3(4c2 + 8C3 + 4c4 + Tcg + 4C7) 

where we have repackaged the terms independent of (x, y, z) in Eq. (22) into the coefficients 
"Ti, T 2 , r 3 by re-expressing 1 as ti -|- r 2 -|- r 3 . Thus, the necessary and sufhcient set of bounds on 
quartic curvature corrections in D = b are 

— 8 C 3 -|- 4c4 cq 2cj > 0 

32ci -|- 8 c2 -|- 8C3 -|- 4c4 -|- 5c6 -|- 2c7 > 0 (35) 

4c 2 -|- 8C3 d- 4 c4 7cq 4c7 > 0, 

coming from analyticity of the four-point graviton scattering amplitude. 
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3.3 Theories in D > 6 


Consider finally the general case of -D > 6. It is a non-trivial linear algebra problem to determine 
the parameter space of (x, y, z) corresponding to physical polarization conhgurations. Each 
physically allowed point {x,y,z) yields a positivity bound via Eq. (22). The set of all positive 
linear combinations of such bounds is given by plugging into Eq. (22) the set of all points in 
the convex hull S spanning physically allowed values of {x,y,z). Fully characterizing all such 
(x, y, z) is beyond the scope of the present work. However, we can derive a general collection 
of necessary conditions from a subset of extremal vertices on the boundary of S. The details of 
the calculation are given in App. A, but the vertices are 


vi = (0,0,1) 

V2 = (l) 1 ~ D-2 ri-3’ 


^3 = ( 0 , 2 ^> 0 ) 


^4 (^1; £)_2 

V5 = ( 0 , 0 , 0 ) . 


1 + 


4{D mod 2) 


,0 


(36) 


These vectors can be realized by physical polarization choices. The bounds associated with 
the (x, y, z) values in Eq. (36) are necessary for analyticity of four-point scattering amplitudes 
and moreover are a subset of the minimal basis of sufficient bounds. Numerical evaluation, via 
the explicit computation of x, y, and .2 for pseudorandom, traceless, unit-norm matrix pairs of 
various dimensions, shows that the convex hull dehned by the vertices in Eq. (36) is in fact equal 
to the full hull S for even D but is slightly smaller than S in odd D. Note that the vectors (36) 
are a generalization of those we saw in earlier sections, so vi, V 2 , and vs coincide with the 
vectors from D = 5. Moreover, each corner corresponds to a certain extreme conhguration of 
polarizations. For example, Vi corresponds anticommuting polarizations as in Eq. (32), while 
V 2 i Vs, V 4 , and correspond to commuting polarizations. For the latter, the polarizations are 
mutually diagonalizable and can without loss of generality be represented as traceless diagonal 
matrices. See App. A for details. 

Plugging the vectors in Eq. (36) into the bound in Eq. (4), we obtain the positivity bounds 


2(1 


3 

D-2 


+ 


1 

D-3 


) 


(sh) 


1 + 


4{D mod 2) 
iD-l){D-3) 


— 8C3 -|- 4c4 Cq 2c7 > 0 

( 4 c 2 -|- 8 C 3 -|- 4 c 4 Cq Cy) -|- 32Ci -|- 4 c 2 + 4cg -|- Cy > 0 

( ^_ 2 ^ ( 4 c 2 -|- 8 C 3 -|- 4 c 4 -|- Cg -|- Cy) -|- 2 c 6 -|- Cy > 0 

( 4 c 2 -|- 8 C 3 -|- 4 c 4 -|- Cg -f- Cy) -|- 32ci -|- 4 c 2 + 4cg -|- Cy > 0 


2cg -|- Cy > 0, 


(37) 
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which are a stringent set of requirements on quartic curvature corrections to general relativity 
in D > 6, necessary to guarantee analyticity of scattering amplitudes. 

3.4 Supersymmetric Theories 

We now consider supersymmetric quartic curvature corrections. Conveniently, Ref. [17] derived 
a basis for independent off-shell supersymmetric quartic curvature invariants, 

A = AOa + BOb + COc, (38) 

where Oa, Ob, and Oc are proportional to more familiar looking forms denoted in the literature 
[17,18] by tstsR'^, and eioeioR^, respectively. In terms of the basis defined in Eq. (4), 

these supersymmetric operators are 

Oa = Oi- 16C>2 + 2C>3 - 32C>5 + WOq + 320r 

Ob = -Oi + 8 C >2 - 2C>3 + 4C>4 (39) 

Oc = Oi- 16C>2 + 2C>3 + I 6 O 4 - 32C>5 + 16C>6 - 32C>7, 
corresponding to the following choice of operator coefficients: 

ci = A-B + C C 2 = -164 + SB - 16^ cg = 24 - 2R + 2C 

C 4 = 45 + 16C C 5 = -324 - 32C Cg = 164 + 16C (40) 

C 7 = 324 - 32 a 

Plugging this choice into Eq. (22), we obtain 

A + B{y + z)>0. (41) 

Note that C drops out of the calculation completely, since at quartic order in graviton pertur¬ 
bations it is a total derivative in all dimensions [13]. Recall from App. A that while simple 
matrix identities imply that y + z < 1, no point in the hull S actually saturates this bound. 
For example, in 5 = 4, Eq. (28) implies that y + z = ^, so A + ^B > 0. In D = 5, Eq. (32) 
implies that \ <y + z<]^, so A + |5 > 0 and 4 -|- |5 > 0. Finally, in 5 > 6 , inputting 
the vertices in Eq. (36) into Eq. (41) yields the complete set of positivity bounds for quartic 
curvature operators in supergravity theories. In summary, we find 

4 + |5>0 (A* = 5) 

4 > 0 (5 > 6 ) (42) 

+ (1 “ ^ > 0) (any D) 

noting that in 5 = 4 and 5 = 5, the final bound reduces to 4 -|- ^5 > 0. 


13 




3.5 String Theories 


As a consistency check, we now apply our bounds to string theory, which is arguably the leading 
candidate for the ultraviolet completion of gravity. Conveniently, quartic curvature corrections 
have been been dutifully computed at tree-level in the existing literature for the bosonic [19,20], 
type II [21-23], and heterotic string [22,23]. The type I string is dual to the heterotic string and 
has the same low-energy effective action [18,24], so we need not consider it as a separate case. 
The resulting effective theory is described by 

A = AO A + BOb + COc + AOa, (43) 

where Oa, Ob, and Oc are the supersymmetric operators from the previous section and Oa is 
a non-supersymmetric operator dehned as 

Oa = —Oi -j- IOO 2 -l- O 4 . (44) 

In various string theories, the operator coefficients are 



A 

B 

c 

A 

bosonic 

C(3) 

0 

-C(3) 

16 

type II 

C(3) 

0 

-C(3) 

0 

heterotic 

C(3) 

1 

-C(3) 

0 


where each entry is normalized by a factor of a'^/1024fi:^. 

As expected, since the type II and heterotic string theories are supersymmetric, their co¬ 
efficients in Eq. (45) satisfy the bound for supersymmetric theories in Eq. (42). Since the 
bosonic string is non-supersymmetric, the bound is more complicated. In particular, plugging 
the corresponding operator coefficients into Eq. (22), we obtain 

C{3) + 2{x + lly + z)>0, (46) 

which is indeed positive, as x,y, z > 0. Thus, we have verihed that quartic curvature corrections 
in bosonic, type II, and heterotic string theory are consistent with unitarity and analyticity. 

4 Bounds on Quadratic Curvature Corrections 

Next, we consider analyticity constraints on C 2 , which characterizes quadratic curvature correc¬ 
tions in the graviton effective theory. As shown in Ref. [9], the Gauss-Bonnet term 

C 2 = + R^) (47) 
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does not introduce ghost modes in any dimension D, so in this basis the graviton propagator is 
unmodihed. To avoid ghost pathologies, we only consider curvature invariants of this form. For 
D = A, the Gauss-Bonnet term is furthermore a total derivative and thus does not affect local 
dynamics. As recently shown [25], however, the Gauss-Bonnet term is critical for computing 
and interpreting the leading ultraviolet divergences of pure gravity. 

Expanding to leading order in the Gauss-Bonnet coefficient A, we compute the quadratic 
curvature correction to the graviton scattering amplitude in the forward limit. 


M2{s,t ^ 0) = AXk'^s'^ 


l-Ll^ fju I fju I fjiU fju 

^3up^2 ^2ap. 




flU 




fll/ 




- [k2k'^€2^€4p^eia^e'^ + ^l^3ei/e3p^e2Q/3G 


a/3 


(48) 


where we have expanded formally in f-dependence arising from propagator denominators, but 
we have yet to evaluate the numerators in the forward limit. 

The hrst line of Eq. (48) is manifestly regular in the forward limit t = 0, so for these terms 
we can simply set €3 = ei and €4 = € 2 - On the other hand, the second line of Eq. (48) is naively 
singular since 1/t diverges as f ^ 0. However, this singularity is canceled by the numerator 
factor, which vanishes in the forward limit as €3 —)■ ei and €4 —)■ € 2 - It will be convenient to 
rewrite this expression in terms of the momentum transfer. 


q = ki + k3 = -{k 2 ki), 


(49) 


where t = —q^. For real kinematics, q is spacelike and vanishes in the forward limit. Note that 
qpq'^/q^ is simply a projection operator in the direction of the spacelike exchanged momentum. 

We note that k^ is simply a real spatial rotation of — fci, and likewise for k^ and k 2 - By 
symmetry, this then implies that e^'^k^y = e^’^kip = q^ and = e^k 2 p. = —at 

leading order in q. Rewriting Eq. (48) in terms of g, we then have 


M 2 (s, f —)■ 0 ) = 4Ak^s^ 


fjjiy pu I r\ po 

^2pa^2 + ^3vp^2 ^2crp 


2qPq’' 


P ^/3 I p Ol^ 

2p, ^2pu^la^^\ \ ^lpu^2a/3^2 


(50) 


which is regular because the projection operator q^q’^/q'^ is hnite in the forward limit. To obtain 
a bound on A, we consider all possible choices for the external momenta and polarizations and 
impose positivity bounds on the forward amplitude in Eq. (50). 

As expected, quadratic curvature corrections to graviton scattering scale as M 2 AkS^ so to 
extract an analyticity bound we should apply Eq. (16) for a second-order residue, corresponding 


15 







to n = 2. Unfortunately, this choice also extracts the t-channel singular contribution from 
leading-order graviton exchange, Mi ~ —n^s^/t. In the forward limit, this contribution is 
formally infinite. Of course, in any physical experiment there is an infrared scale p that regulates 
these contributions from long distance physics. This would arise, e.g., from a finite detector 
resolution or beam width [26]. As is common practice for infrared divergences in scattering 
amplitudes, we introduce a mass regulator, sending t ^ f in the denominator. This 

approach was also used in Ref. [1] to make sense of a theory of interacting massless scalars with 
trilinear couplings. Note that as in gauge theory, the mass /i^ is a formal regulator that leaves 
the number degrees of freedom untouched—so the vDVZ discontinuity [27,28], which arises for 
a physical graviton mass included via a Fierz-Pauli Lagrangian term, does not apply here. 

While tames the formal infrared divergences, for A < 1 the forward amplitude will be 
dominated by finite but large contributions from Einstein-Hilbert interactions because |Mi| ^ 

I M 2 1 in this regime.® However, by explicit calculation, we can see from Eq. (18) that Mi ~ 
which is positive. So while positivity is satished, we learn nothing beyond what is 
already borne out from scattering via the leading Einstein-Hilbert term. 

To place a bound on the coefficient A, we must then restrict to a parameter regime where 
I Ml I < I M 2 1, so the contributions from graviton exchange are subdominant to those from the 
Gauss-Bonnet term. This implies that l//i^ < |Afi:^|. Together with the requirement that 
|s| 3> necessary to treat p as a regulator, this forces us to consider the regime 

> p > A, (51) 

where A ~ is the scale of the would-be natural cutoff associated with the derivative 

2 

expansion. We assume throughout that A go that it is below the Planck scale in D 

dimensions. 

Of course, Eq. (51) points to a naively pathological region of the effective field theory, given 
the reasonable expectation of new degrees of freedom of mass m where m ~ A. Moreover, 
Eq. (51) indicates that the infrared regulator p must be larger than some other energy scale 
A. Nevertheless, one can a priori envision an ultraviolet completion in which m 3> A, so new 
degrees of freedom enter at a parametrically higher scale. In that case, A is not the scale of 
any physical states in the theory and is merely the combination of parameters that appears in 
the higher-dimension operator. Indeed, at the level of the scattering amplitude, there are no 
discontinuities that appear around A to signal new degrees of freedom. 

^Note that taking t strictly to zero is not required to derive a positivity bound [29] and positivity holds for any 
non-negative t below gP' [30]. However, we will not need this more general result for our purposes. 


16 



Thus, /i remains smaller than any physical mass scale in the theory and indeed can be treated 
consistently as an infrared regulator. In the absence of new states at A, the Gauss-Bonnet term 
acts effectively as a primordial contact operator over a wide range of scales. Precisely such a 
scenario was considered in Ref. [15], where it was found that such a low-energy effective theory 
is acausal without new states at or below m < A. Other authors [31] have likewise argued that 
a pure Gauss-Bonnet theory is inconsistent with black hole thermodynamics. We will likewise 
hnd a pathology in this theory coming from unitarity and analyticity. 

To apply constraints from unitarity and analyticity, we must hrst ensure that the low-energy 
theory is sensible enough that we can even speak of a long-distance scattering amplitude. Indeed, 
Eq. (51) is plainly strange since |s| 3> A^ violates the derivative expansion. This was required in 
order for the Gauss-Bonnet interactions to dominate over the Einstein-Hilbert action, as was also 
assumed in Ref. [15]. Naively, one would expect a gross departure from perturbative unitarity, 
e.g., probability amplitudes greater than one as well as a breakdown of the loop expansion. 
Nevertheless, there is a wide range of scales where neither sickness actually arises. This hinges 
on the fact that the theory depends on A as well as k, the gravitational coupling constant. 

In particular, note that amplitudes can still be perturbatively small in the regime specihed 
by Eq. (51). For example, M 2 ~ M is still sensible provided k is sufficiently small, corre¬ 
sponding to the weak gravity limit. We can make this more precise by considering the leading 
effect of the Gauss-Bonnet term, which is a cubic vertex of the schematic form Insert¬ 

ing this vertex into low-energy amplitudes, we hnd that the theory remains under perturbative 
control provided Xk^ (~ «:/A^) times the appropriate powers of energy is sufficiently small. In 
D dimensions this implies that 




(t)'" 


(52) 


to safely reside within the regime of perturbativity.® Moreover, Eq. (52) also ensures a pertur¬ 
bative loop expansion, since radiative corrections always introduce additional insertions of the 
Gauss-Bonnet interactions. 

For our purposes, we assume a weak gravity limit dehned by Eq. (52), so the low-energy 
theory is perturbatively unitary. When then apply the method of Sec. 2, where the contour 
around the origin in the complex s plane is widened so as to satisfy Eq. (51), ensuring that s 
is large compared to the scale of the infrared cutoff and that the Gauss-Bonnet term dominates 
the amplitude. Note also that the initial contour encircles a region below the heavy particle 


similar statement applies to pions, which have quartic vertices of the form /XXv^ where v is the breaking 

scale and A controls the derivative expansion. The theory maintains perturbative control provided s <C Xv. 
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threshold, m S> A. 

To see how a pathology arises, it will be convenient to dehne coordinates transverse to the 
incoming momenta, (xi, X 2 ,.. .xd- 2 )- Withont loss of generality, we take the forward limit such 
that the infinitesimal momentum transfer lies in the xi direction, which we henceforth refer to 
as the “direction of approach.” In turn, q^q'^/q^ is a projection operator onto this direction. 
Next, we dehne a particular subset of polarizations in the transverse plane, dehned by rank-two 
diagonal matrices of the form 

Xi Xj 

= ^diag( 0 ,..., 0 ,'^, 0 ,..., 0 ,^, 0 ,...), 

with zero entries except in the x* and xj directions. As the only preferred direction is xi, labeling 
the direction from which we approach the forward limit, the relevant physical polarizations are 
g^( 2 , 3 ), The forward limit of the quadratic correction to the graviton scattering 

amplitude in Eq. (50) for various polarization combinations is 

r 0 , 61 = 62 = 

M{s,t^0)= 2AkV X I 4, ei = e 2 = d(2,3) ( 54 ) 

[ — 1 , 61 = and 62 = 

In the hrst case, ei = 62 = corresponding to polarizations that have support in the 

direction of approach. In the case of iA = 4, this is required because the transverse space only 
has two dimensions. As expected, the amplitude vanishes in this regime because the Gauss- 
Bonnet term is a total derivative in iA = 4. Meanwhile, the second case, ei = 62 = occurs 
when both polarizations are orthogonal to the direction of approach. Of course, this requires 
dimensions D > 5. Finally, in the last case, ei = and 62 = the polarizations occupy 

different planes but share support in the direction of approach, which is only possible in D > 5. 

The upshot of Eq. (54) is that in iA > 5, different polarization conhgurations can yield 
opposite signs for the corrections to the forward scattering amplitude. As a result, this excludes 
both signs of A and thus forbids it entirely. Of course, we made the assumption of Ref. [15] that 
the Gauss-Bonnet term is an effectively primordial contact operator insofar as new degrees of 
freedom enter only at a scale far above the naive cutoff. Hence, the positivity violation in Eq. (54) 
simply implies that this assumption is false. We conclude that a primordial Gauss-Bonnet term 
is forbidden and new degrees of freedom are required at or below the cutoff A. 
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5 Conclusions 


In this paper, we have derived rigorous bounds on the coefficients of quartic and quadratic 
curvature corrections in the low-energy effective theory of gravitons. Our results hinge on 
very general principles: quantum mechanical unitarity and analyticity of scattering amplitudes. 
Consequently, these constraints apply to any consistent perturbative ultraviolet completion of 
gravity. For the quartic curvature operators dehned in Eqs. (3), (4), (24), and (25), we derived 
the positivity bounds in Eq. (30) in D = 4, Eq. (35) in D = 5, and Eq. (37) in arbitrary 
D > 6. We also presented constraints on supergravity theories and checked that all of our results 
are consistent with known calculations in weakly-coupled string theories. For the quadratic 
curvature correction in Eq. (47), we showed that both signs of its coefficient A are inconsistent 
unless new degrees of freedom enter at the natural cutoff A ~ specihed by the effective 

theory. In short, a primordial Gauss-Bonnet term is forbidden. 

Many possibilities remain for future work. While four-point graviton scattering cannot probe 
curvature operators beyond quartic order, little is known of higher-point amplitudes. Such ampli¬ 
tudes are functions of many more kinematic invariants and should thus enforce commensurately 
more positivity constraints. Another issue meriting further study is that of cubic curvature 
operators. Here, positivity bounds encounter technical challenges due to the vanishing of the 
associated amplitude in the forward limit [30,32]. As noted in Ref. [32], this problem is closely 
related to the a-theorem in iA = 6. 

Distinguishing low-energy effective theories that are consistent with ultraviolet completion 
from those that are not presents a signihcant challenge. Systematizing this procedure is im¬ 
portant for delineating the space of possible physical laws, but has also become important for 
model-building in more phenomenological contexts [33] and in inflation [34-38]. In this paper, 
the low-energy tools of analyticity and unitarity enabled us to hnd solutions to this problem 
in gravitational theories, allowing us to constrain higher-curvature corrections to gravity in our 
own universe—applying our quartic curvature results to D = 4—and further discover bounds 
applicable in any consistent theory. 
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A Bounding Invariants in General Dimension 

We have shown that the graviton scattering amplitude can be expressed in terms of invariant 
products of graviton polarizations ei and € 2 , which are real, symmetric, traceless {D — 2)-by- 
{D — 2) matrices with unit normalization. To recapitulate from Sec. 3, given the Hermitian 
matrices 


H+ = {ei,e 2}/2 and H_ = i[ei,e 2 ]/ 2 , (A.l) 

we can define the invariants 


X = Tt{H+)Tt{H+), y = TT{H+-H+), z = H_). (A.2) 


The space of physical polarizations ei and 62 then maps onto a physical region in {x, y, z), which 
through Eq. (22) implies positivity constraints on operator coefficients in the effective theory. 

What are the bounds on {x, y, z)7 We first note that since and are Hermitian, their 
squares are positive semidehnite, so x,y,z > 0. Moreover, a straightforward application of the 
Cauchy-Schwarz inequality implies x = Tr(ei • e 2 )Tr(ei • 62 ) < Tr(ei • ei)Tr(e 2 • € 2 ) = 1, with 
equality if and only if ci = ± 62 , and similarly y + z = Tr(ei • Ci • 62 • 62 ) < Tr(ei • ei)Tr(e 2 • 62 ) = 1. 
There are, however, many additional constraints on {x,y,z), which we now discuss. 

Crucially, a weighted average of any number of positivity bounds yields another valid posi¬ 
tivity bound. This implies that a space of necessary conditions can be derived by constructing a 
convex hull S in {x, y, z) that contains the physically allowed region. Without loss of generality, 
S is 


I/, ^) = < X] 


TiVi 


i=l 


Tj > 0 and Ti = 1 


(A.3) 


i=l 


where n* denote extremal points. In this Appendix, we will construct the subset of the Vi that 
are on the edges of the unit cube in {x, y, z); let the convex hull described by these vertices be S. 
In even dimension, numerical evaluation suggests that S = S, while in odd D > 6, it is possible 
for points to he slightly outside S. 
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Let us first consider the case where ei and €2 are anticommuting, so x = ?/ = 0 and we wish 
to maximize 2 ;. Going to a basis in which ci is diagonal, we hnd that antisymmetry of ei • €2 
implies that for each i, j, 

+ ^ljj)^2ij = 0 (A.4) 

and 

2 = -Tr(ei ■ € 2 -€i- 62 ) = Y^ 4n4ij- (A.5) 

hj 

Since Eq. (A.4) implies eiiie 2 ii = 0 for each z, the normalization condition 1 implies 

by Eq. (A.5) that nonzero diagonal terms in €2 can only decrease z. We therefore take €2 to have 
vanishing diagonal. Similarly, since Yhi 4ii fixed to unity, we should require that, for each i for 
which eiii 7 ^ 0 , there exists j such that e 2 ij 7 ^ 0 ; letting be nonvanishing for some Zq even if 
^ 2 ioj = 0 for all j would decrease 2 ; by Eq. (A.5). Writing 4ij = Pi^ where J2iPi = 1 by fh® 
normalization constraint, we can then consider z to be a weighted average over the Thus, 
is maximized when we weight the average most in favor of the z for which is maximal. 
Suppose there are N such z, which we can without loss of generality take to be 1 through N, 
for which takes its maximal value, z.e., = max* efjj = for all i* G {1,..., A^}. Then 

is maximized when we have pi = 1/N for z G {1,..., A^} and Pi = 0 otherwise, for which we 
obtain z = . Finally, it remains to determine the maximal possible value of . Since e\ is 

of unit norm, its maximal value is attained when we load as much of the normalization into as 
few of the as possible. By tracelessness of ei, at least two of the em must be nonzero. Thus, 
takes its maximum value of 1/2 when ei oc <73 in some 2 -by -2 block, up to permutation of 
coordinate labels. That is, a choice of polarizations that maximizes 2 ; for x = z/ = 0 is 

ei = ^ 0-3 © On -4 and €2 = © On- 4 , (A. 6 ) 

which yields the point 

^^i = (0,0,i). (A.7) 

Let us henceforth consider the case where z = 0 and explore in x, y. This means that the 
(real, symmetric) matrices 61,2 commute and so are simultaneously diagonalizable. Taking x = 1, 
we can ask how large y can be, which will give a vertex of S. Since Ci = ±62 for x = 1, we 
have y = Tr(ei • ei • ei • ei). That is, y has positive hrst and second derivatives in each of the 
\eiii\ values; y is therefore maximized when one of the \eiii\ is as large as possible and the others 
are equal and small. (If the smaller numbers in the list were unequal, we could always make y 
larger by shifting some weight back to the element in the list with the largest absolute value.) 
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That is, a choice of polarizations maximizing y for x = 1 and z = 0 is 


= 7(5d)(5r!)'“‘‘S(^' 1. 3)). (A.8) 

which corresponds to the vertex 

= (l, 1 - + ;^,0) . (A.9) 

Next, still taking z = 0, we consider a different extreme, setting x = 0 and maximizing 
y. Again simnltaneonsly diagonalizing ei and 62 , we have y = Analogously with 

the case of Ui, we can write pi = and consider |/ to be a weighted average over the e^a- 
Let maxj and, without loss of generality, suppose that e^u = for i G {1,..., A^} for 

some N. Then y is maximized if we take pi = 1/N for i G {1,..., A^} and Pi = 0 otherwise, in 
which case y = Now, by the unit normalization of 62 , is maximized when as much of the 
normalization as possible is loaded into as few terms as possible, he., N is minimized. Since ei 
is traceless, at least two of the pi are nonzero, so N >2. The maximum value of thus occurs 
when N = 2, which fixes ei oc cxs © ^D-i- We now must maximize the common absolute value 
of the hrst two entries in e 2 M, subject to the constraints that Yhi = 0 , ^ 2 ii = 1; since 
a; = 0, ^iii^ 2 ii = 0. This last constraint implies that the first two entries in e 2 ii have the same 
sign. Thus, y is maximized for x = z = 0 for the choice of polarizations 

© Oli -4 and 62 = y^^^z^l^^diag (^, -1,..., -l) , (A.IO) 

for which we hnd the vertex 

’’3=(0.2g©.0)- (Vll) 

Note that the polarization configuration in Eq. (A.IO), and hence the vertex in Eq. (A.11), 
requires D > 5. 

On the other hand, we can minimize y/xfoi z = 0. Using symmetry and reality to diagonalize 
7L+ = diagh+, we have x = |h+ • np, where n = (1,1,..., 1), so the Cauchy-Schwarz inequality 
implies 

X < |h+|^|n|^ = {D — 2)y. (A.12) 

Eq. (A.12) is saturated when oc Id- 2 - If D is even, this choice is possible with 

*^1 = ^2 = ^T^diag (1, -1,1, -1,...), (A.13) 

which yields 

{x,y,z) = . (A.14) 
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Let us now consider the odd-dimensional case where z = 0 and x = 1. Simultaneously diag¬ 
onalizing ei and € 2 , we have y = Again, y has positive first and second derivatives in 

|eiij|, so it is minimized when the |eijj| are all equal. In odd dimension, this is not possible 
while retaining tracelessness, so the best one can do, making the |eiji| as equal as possible, is 
the choice 

£1 = £2 = • — ; \ i •••) (A. 15 ) 

D-1 

2 

which results in the vertex 


{x,y,z) ( 1 , + £)_3 D- 2 ’^)- 


(A.16) 


We can combine Eqs. (A. 14) and (A.16) to write the vertex of S as 


U4 


f 1 1 -I 4(D mod 2) 

D-2 {D-l){D-3) 



(A.17) 


We note that for both D = A and D = 5, V 2 and U 4 are the same point. Moreover, 
generalizes the third vertex applicable in H = 5, while the polarization choice for U 3 does not 
apply in H = 4. In H > 6 , there is one remaining linearly independent vertex, which can be 
obtained by choosing ei ■ €2 = 0 ^- 2 , e.g'.. 


^ ;^(T 3 © 0 d _4 and 62 ^ ffi <^ 3 , 


(A.18) 


which results in the point 


!l5 = (0,0,0). 


(A.19) 


Together, Eqs. (A.7), (A.9), (A. 11 ), (A.17), and (A.19) are the vertices of S given in Eq. (36). 
They correspond via Eq. (22) to a set of linearly independent bounds (37) that must be satisfied 
in any gravity theory in H > 6 . 
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